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Abstract 

In the article correct method for the kinetic Boltzmann equation asymptotic solution is formu- 
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method, within Enskog's approach, are equivalent up to the infinitesimal first order terms of the 
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lent gas flow is proposed, as stratified on components gas flow, which is described by the derived 
equations system of multicomponent non-equilibrium gas dynamics. 
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I. INTRODUCTION 



In 1912 Hilbert considered the kinetic Boltzmann equation for one-component gas as 



an example o: 
solution (see 



' integral equation and proposed a "recipe" for its approximate (asymptotic) 
l), Chapter XXII). Hilbert 's "recipe" was inconvenient for practical use, be- 
cause the five arbitrary functional parameters of the first and the following approximations 
of the velocity distribution function had to be found by solving the differential equations 
in partial derivatives (equations of gas dynamics of the first and higher orders). Five years 
later Enskog proposed to use zero conditions, conditions (I57p - (|B^|) below with zero right- 
hand sides, to determine the five arbitrary functional parameters of the first and following 
approximations of the velocity distribution function. As we will see in the section [Til (in case 
of one-component gas), the imposition of the zero conditions leads, in fact, to using different 
comparison scales in the asymptotic expansion of the velocity distribution function and in 
the asymptotic expansion of the particle number density, the mean (mass) velocity and the 
temperature, that are derived from the asymptotic expansion of the velocity distribution 
function by integration over velocities with different weighting functions. As a result of par- 
alogism of the method of successive approximations (one has to set variable coefficients of 
the same terms of the unified comparison scale equal to each other) partial time derivatives 
vanish in the necessary conditions of solutions existence of integral equations of higher or- 
ders (see below) and with them terms of gas-dynamic equations, corresponding to viscosity, 
heat conduction, . . . vanish. Enskog "improved" the situation by the introducing (see, for 
example, [2], Chapter 7, § 1, Section 5) of the unsubstantiated expansion of partial time 
derivative. 

The approach of Struminskii, who had proposed in 1974 in |3] his approximate (asymptotic) 
method of solution of the system of kinetic Boltzmann equations for multicomponent gas, 
differs from the approach of Enskog to asymptotic solution of the Boltzmann equations sys- 
tem for gas mixture in that, how the infinitesimal parameter is introduced in the Boltzmann 
equations system for gas mixture, i.e. the solution is constructing in another asymptotic 
limit. It could be noted, that approaches, that are similar to Struminskii's approach to 
asymptotic solution of the kinetic Boltzmann equations system, were regarded earlier in the 



kinetic theory of plasma (for instance, see j^], § 7.5). In substance, Struminskii's method of 
solution of kinetic equations system is the same as Enskog's method (Struminskii used the 



partial time derivative expansion, as Enskog did). 

After some preliminary remarks in section [Til m section UTTI there will be proposed the cor- 
rect method of asymptotic solution of the kinetic Boltzmann equations system for multicom- 
ponent gas mixture for the approach, that combines Enskog's and Struminskii's approaches; 
in particular, it will be shown, how one has to modify Enskog's method - in addition to 
asymptotic expansion of the velocity distribution function of z-component particles of gas 
mixture it is necessary to determine and to use the expansion of the particle number density 
Hi of i-component, mean mass velocity u and temperature T of the gas mixture. 

Further, in the section [TV] the system of infinitesimal first order equations of multicom- 
ponent n on- equilibrium gas dynamics, appearing during the process of the solution of the 
system of Boltzmann equations by successive approximations method in the section II III as 
necessary condition of the existence of approximate (asymptotic) solution of the integral 
equations system, is considered in more detail. 

In sect ion \V\the expressions for the velocity distribution function of particles of gas mixture 
in the infinitesimal first order are obtained for the approach, that Enskog's approach and 
Struminskii's approach combines, in particular it is shown, that the velocity distribution 
functions, that are derived within Enskog's approach by Enskog's method and the proposed 
method, are equivalent up to the infinitesimal first order terms of the asymptotic expansion 
(inclusive), but, generally speaking, differ in the next order of the asymptotic expansion; 
the difference is a possible reason, why the transition to a higher order approximation in 
Enskog's method does not lead to any essential improvements of the result. 

In section I VI I an interpretation of turbulent flows is proposed within the multicomponent 
gas dynamics. 

Struminskii's article 3], and, as a result, the followings articles, that referred to the 
paper jsj (for example, p]), contained errors in calculation of collision integrals. In appendix 
IA1 general analytic expressions of these collision integrals (that depend on the interaction 
cross-section) are derived in general case, when separate components (with Maxwell velocity 
distribution function of particles) have different mean velocities and temperatures. 

For the elementary model of an interaction of gas mixture particles - the interaction 
of particles as rigid spheres - collision integrals can be simply calculated completely, the 
expressions for the interaction potential of rigid spheres are written in appendix Q5] 

Notations, used below, are close to notations in Q; it is assumed, that all regarded 



functions are continuous and continuously differentiable so many times as it is necessary, if 
their derivatives are considered, and all regarded integrals converge. 



II. SOME REMARKS 

In kinetic theory of gases and gas dynamics, closely connected with it, there are no small 
dimensionless constants, similar, for example, to the fine structure constant in quantum 
electrodynamics, in terms of which it would be possible to construct perturbation theory 
expansions, paying special attention to convergence of derivable expansions. In particular 
applications of kinetic theory of gases and gas dynamics it is possible often to make up from 
parameters of the problem a variable dimensionless value (in gas dynamics such variable 
dimensionless values are named numbers: the Mach number, the Reynolds number . . . ), we 
will denote this value by 9, and to consider asymptotic expansions with variable coefficients 
([6|, Chapter V, § 2, Section 5) of physical values in terms of integral powers of 9 - {9 r } r>0 , 
when 9 tends to (most important asymptotic expansions, further just those asymptotic 
expansions are considered; generally speaking, it is possible to use also integer negative 
powers r, but we shall not do it), or more rarely {9~ r } r>0 , when 9 tends to +oo or to — oo; 
changing 9 to 1/9, it is possible to switch easily from one of these cases to another. In 
the theory of asymptotic expansions variable coefficients are considered to be comparable 
to the unit (and may depend on 9 in general case, but not in this article) and functions of 
a comparison scale {{9 r } r>0 ) - to be infinitesimal. If an asymptotic solution to specified 
orecision (9 r ) relative to a comparison scale ({9 r } r>0 ) exists, it is unique (see, for example, 
6|, Chapter V, § 2, Section 2) and, if asymptotic expansions are considered at 9 — > 0, is 
exact, at least, at 9 = 0, therefore it is convenient to consider the asymptotic expansions at 
9 — > (and integer r > 0). 

In applications of the theory of asymptotic expansions, solving physical equations by 
the successive approximations method, it is convenient not to extract the dimensionless 
parameter 9 r from a term of asymptotic expansion X^j cLS cL really small multiplier, but 

to use 9 r as an indicator of infinitesimal order of corresponding (infinitely) small variable 
coefficient (x), setting in the resultant expressions 9 = 1; with such agreement the 
formal "small" parameter 9 can be introduced in an equation arbitrarily, but the meaning 
of obtained asymptotic solutions of such equation with small parameter is determined by 



4 



physical validity of the introduction of the small parameter in the original equation. 

By method of successive approximations we obtain, generally speaking, an asymptotic 
solution of the problem. If some additional conditions exist (for example, initial conditions 
or boundary conditions for a differential equation), the conditions must be expanded in 
asymptotic series, and the obtained equations system have to be solved, setting variable 
coefficient at the same terms of unified comparison scale equal to each other. If the system 
is solved successfully, we obtain an asymptotic solution of the initial equation, that meets 
(asymptotically) the additional conditions. In some cases the found asymptotic solution 
turns out to be regular solution (Q, Chapter 1), i.e. it depends analytically on 6, or even 
exact solution of the problem. For example, the equation (see js|, Chapter V, § 2 or |9|, 
Chapter IV, § 7.1) 

#f + / = 0, /(0) = (1) 

with zero initial condition has asymptotic solution (t) =0, (r = 0, 1, 2 . . .) relative to 
the comparison scale, formed of functions 6 r , that is also the exact solution of the problem. 
But the problem 

#f + / = 0, /(0) = 1 (2) 

has no analogous asymptotic solution, because the zero order asymptotic solution of the 
differential equation (i) = contradicts to the initial condition /(°) (0) = 1. It is not a 
serious drawback of the method of successive approximations. In the problem (T5]) one can 
introduce new function g (t): 

f(t)=exp(-t/9)g(t) (3) 

(cf. with [7]; though not all considered by Lomov expansions are the asymptotic expansions 
with variable coefficients, and, accordingly, Lomov's approach in toto seems not quite correct, 
in the monograph [7] it is actually shown, that similar replacements of functions allow to 
obtain asymptotic solutions for a wide class of problems) and obtain asymptotic solution 
of the problem by the successive approximations method g^ (£) = const = 1, g^> (t) 
const = (r = 1, 2 . . .), that is again the exact solution. 

Criticism of the successive approximations method in [8|, Chapter V, § 2 and [9|, Chap- 
ter IV, § 7.1 reflects, possibly, dissatisfaction of the authors with the introduction of the 



unsubstantiated expansion of partial time derivative in Enskog's method for the kinetic 
Boltzmann equation asymptotic solution, see below. 

Hilbert noted in [1] , Chapter XXII, that the expansion of the velocity distribution function 

F = y + * + + • • • (4) 
A 

[similar to (|34|) below; Hilbert considered one-component gas only; we use Hilbert's notations 
here, however, the meaning of the notations is sufficiently understandable from context] is 
the power series in (a small parameter) A, satisfying to the Boltzmann equation and such, 
that expressions [cf. with (!50l)-( !53l) and (j54~l) - (l56l) below] 

J ip^F du=jj du + J dco 

+ A J \p^Xdu + --- (2 = 1,2,3,4,5) (5) 

for t = to change to power series 

f(i) 

A {1) = ±—+g^+\h^ + --- (2 = 1,2,3,4,5), (6) 
A 

and in the theorem, completing the work, he has formulated a "recipe" to get an asymptotic 
solution for the Boltzmann equation, in which he proposed to determine five arbitrary func- 
tional parameters of functions $, ^f, X . . . " from five partial differential equations 11 [analogous 
fl82|) . (|99D below], "setting at t = t " 

J dcu = AA (i) (2 = 1, 2, 3, 4, 5), (7) 

Jij^^du = (2 = 1,2,3,4,5), (8) 

ip (i) Xduj = (i = 1,2,3,4,5). (9) 

In notations from (HSj) . fl50|) - fl59|) below, Hilbert proposed simply to set special initial values 

n(r,t ) = n^(r,t ), (10) 
u(r,t ) = u^(r,t ), (11) 
T{r,t ) = T®{r,t ) (12) 



or 



i) {l) f {r) dc 



= (/ = 1,2,3) (13) 

t=to 
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for r = 1 , 2 . . . . 

"For the further substantiation of the gas theory" it might be as well to supplement 
Hilbert's theorem with explicit expressions of five arbitrary functional parameters of the 
functions , found on the r-step (r = 0, 1, 2 . . .) of the successive approximations method, 
through through the physical parameters of the gas [see (158]) - (190]) . (jl05p - (jl07p below]; Hilbert 
had not done it. 

Enskog supplemented Hilbert's "recipe" . However, at that Enskog made a logical mistake. 
He used zero conditions (TT3T) identically at any t, not at t = to only (see [2j, Chapter 7, § 1, 
Section 1): 

yV'¥ r) dc = (1 = 1,2,3) (14) 

for r = 1, 2 . . . . From the viewpoint of the theory of asymptotic expansions, Enskog [instead 
of (|54|) - fl56|) below] supposed 

n (r, t, 9) = 9° n (r, t, 9) + 9 1 + 9 2 + • • ■ , (15) 
u(r,t,9) = 9°u(r,t,9) + 9 1 + 9 2 + -- - , (16) 
T (r, t, 9) = 9° T (r, t, 9) + 6 1 + 6 2 + • • • . (17) 

If n, u and T did not depend on r and t, that would mean, that Enskog used different scales 
of comparison in the method of successive approximations at the same time {n (6) , 6 1 , 8 2 . . .}, 
{u (6) , 9 l , 9 2 . . .}, {T (8) , 9 1 , 9 2 . . .}, that is in itself wrong. In general case, when n, u and 
T depend on r and t, sums (|T5|) - (|T7]) can not even be considered as asymptotic expansions 
with variable coefficients. 

Violation of logic of the successive approximations method shows itself immediately in 
vanishing of partial time derivatives in gas-dynamics equations systems of the (r + l)-order 
[analogous fl82|) (r = 1,2, . . .) below, the necessary conditions of solutions existence of the 
integral higher orders equations in Enskog' s approach], according to (Tbfl) . 

J 4> {l)d -^dc ^)fir) dc = o (i = h 2) 3), (18) 

and with them terms of gas-dynamic equations, corresponding to viscosity, heat conduction, 
. . . , vanish. 

Hence Enskog should have concluded, that the use of (fT4|) is incorrect, as in proof of 
theorem by contradiction, but instead of that Enskog proposed to use the unsubstantiated 
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expansion of partial time derivative 



III. CORRECT METHOD OF SOLUTION OF THE KINETIC BOLTZMANN 
EQUATIONS SYSTEM 

The Boltzmann equations system, that describes change of dependent on t and spatial 
coordinates, prescribed by radius-vector r, the velocity distribution functions /j (t, r, Cj) due 
to collision with particles of other components of mixture of rarefied monatomic gases, where 
Cj are the velocities of particles of z-component of the mixture {see [2], Chapter 8, (1.1); 
discussion of the derivation of the Boltzmann equations system and its applicability range 



see, for example, in [2|, Chapters 3 and 18, 10], Chapter 7, § 1, |ll|, Chapter 3 and the 



Bogolyubov paper [12] ; below the central interaction of molecules are considered only, when 
the force, with which each molecule acts on the other, is directed along the line, connecting 
the centers of the molecules}, could be written as: 



t+C'-f+^/ir^lll {f - f '- ms -> bdbdedc 

1 1 j€N 




= E // Wi - m) k * dkdc > G N ^ ( 2 °) 

in (120]) iV is a set of indexes, that are numbering components of the mixture; Xj is an external 
force, which acts on the molecule of the z-component; rrii is the mass of the molecule of the 
z-component; is the modulus of the relative velocity of colliding particles = Cj — Cj\ b 
is the impact distance, e is the azimuth angle, k is the unit vector, directed to the center of 
mass of the colliding particles from the point of closest approach - see Chapter 3, Fig. 
3; the scalar function kij (fl^k) is determined by equality 

gij b alb ale == kij (ik; (21) 

by prime in (120]) and below the velocities and the functions of velocities after the collision 
are denoted. 

Let us introduce following notations: 

Ji (fi,f)= ft (fif-mkdkdc, (22) 



Jid (Ji, fj) = jj (fifj - flfj) kij dk dr. - (23) 



to differ velocities of colliding molecules of the same kind in ( 122]) the one velocity is denoted 
by Cj and the other is denoted by c (without any index) and the index of the corresponding 
velocity distribution function / is omitted. 

In Enskog's approach the differential parts of the Boltzmann equations (l20p . that are 
denoted by Pj/j below, are considered to be small as compared with the right-hand sides 
of equations (120]) - see [2j, Chapter 7, § 1, Section 5 - therefore the indicator of infinite 
smallness 9 is formally introduced in the Boltzmann equations system in the following way: 

9V i f i = -J2JiAfiJj) (*eiV). (24) 
j 

In Struminskii's approach to the asymptotic solution of the Boltzmann equations system the 
differential parts of the Boltzmann equations (120]) and the collision integrals of the particles 
of i-component with the particles of the other components are considered to be small as 
compared with the collision integral of the particles of i-component between each other, 
therefore the indicator of infinitesimality 9 is introduced in the Boltzmann equations system 
in another way: 

mfi = ~Ji (fi, f)~ 9 Yl (/" & (« G N). (25) 

It is possible to combine Enskog's approach with Struminskii's approach. For this purpose 
we divide the set of mixture components N into two subsets: the subset of components, that 
we call formally inner components (we could consider the case, when there are some subsets 
of inner components, but this case does not fundamentally differ from the one, considered 
below, the only difference is that the notation become more complicated) and the subset 
of components, that we call external components. To differ the two groups of mixture 
components we denote the subset of indexes of inner components N as well as the indexes of 
inner components i G N by symbol "~" and the subset of indexes of external components N 
as well as the indexes of external components i G N by symbol " ~" ; the intersection of the 
sets N and N is the empty set - N (1 N = and the union of these sets is the set of indexes 
of all mixture components N U N = N; if an assertion concerns both kinds of components 
the symbols """ and """ will be omitted. In new notations the Boltzmann equations system 
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can be rewritten as: 

Mtft = -Y,JQ(f{,fi) (< e N) , (26a) 

9V.fi = ~Ji (fh f) ~ J2 (I,, h) (i 6 N) . (26b) 

Let us write the asymptotic expansion of the velocity distribution function fa of particles 
of i-component as formal series of successive approximations in powers of 9: 

fi = fl 0) + + Pf® + ■■■. (27) 

The differential parts of the equations f[2"6"j) are written as: 

= Vf ] + 9vf + 2 pf ) + • • • , (28) 



where 

A W = ^ + C ,-^ + --^ (r = 0,1,2...), (29) 
at or mi oci 

- cf. with (2!, Chapter 7, § 1, Sections 4, 5 and In (r28|) -( l29|) the partial time derivative 
expansion (TTO]) is not used in contrast to that, how it was made by Enskog and further 
by Struminskii. As result, described below method for solution of the system of kinetic 
Boltzmann equations differ fundamentally from Enskog's method and Struminskii's method. 

Substituting ( 127|) and ( 128]) in (|26aj) and equating coefficients at the same powers of 9 to 
each other, we obtain the equations system of the method of successive approximations for 
finding the velocity distribution functions of inner components particles of gas mixture f? r ^; 
taking introduced notations f[2"2"j) . f[2"3"j) and (j29"|) into account, the system can be rewritten 
as: 

£^(4 0) >if) = ° ( /G 4 (30) 
*>r> + E Jy (4 r) ./f ) +EE% (r s> . + E % (4 0) > /f ) 

+ £ £ {ft*'*' ff) = (« G r = 1, 2 . . .) . (31) 

jev s =o 
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Similarly substituting ( p7|) and ( 128|) in ( 126bl) and equating coefficients at the same powers 
of 9 to each other, we obtain the equations system of the method of successive approximations 
for finding the velocity distribution functions of particles of external components of gas 
mixture ft . 

^(/r,/ (0) )=o (ietf), (32) 

s=l 

+ EE J ( ? (^" 1 ~ S) '^ ) ) =0 (^iV, r = l,2...). (33) 

j^' s=0 

It could be noted also, that in original Hilbert's paper jl| and Enskog's paper (see. 
Chapter 7, § 1, Sections 4) the formal parameter 9 was not directly introduced in 
the Boltzmann equation, but in the series of successive approximations for the velocity 
distribution function the parameter 9 was introduced in another way (not as in (l27p): 

if' E = ^ = ^if + /fW ) + ---; (34) 

Obviously, with the same result it is possible to use the expansion ( 1271) for the velocity 
distribution function, but to enter the multiplier 9 in the Boltzmann equation, as it was made 
in (1241) . Successive approximations f^, f\ l \ fj® . . ., that are calculated within Enskog's 



approach, turn out to be (inversely) ordered by the number density of molecules of the 
mixture n: is proportional to n, does not depend directly on n etc. Thus, there is a 
physical justification to use the method of successive approximations to find an asymptotic 
solution of the Boltzmann equation. For Struminskii's approach it is difficult to determine 
explicitly a small physical variable to construct asymptotic expansions in powers of the 
variable. 

Speaking about an order of approximation below, we assume the order to be equal to 
the value of index r in (1311) . (1331) . According to (1231) . (1301) . in zero order approximation we 
have the following system of integral equations to find the velocity distribution functions of 
particles of inner components of gas mixture fi°^: 

e or. /f ) = e// (/f /f - j? j?) =» a^y (») 

11 



In the kinetic theory of gases the following equality is often used: 

ififjdijbdbdedcidcj = //// ^'ihfj gij b db de dci dcj (36) 




- in ( 13"6"|) . as well as in f ]2"U]) velocities and velocity functions of particles after collision are 
denoted by prime. To each collision of two gas molecules, that transforms the velocities of the 
molecules before the collision Cj, Cj into c/, cj, the collision of two gas molecules corresponds 
bijectively, that transforms the velocities of the molecules c/, cj — >■ — c/, — cj — >■ — q, — Cj — > 
Ci,Cj, and for considered central interactions it follows from the energy conservation law, 
that magnitudes of relative velocities of molecules before and after the collision are equal 

9ij = glj, (37) 

from the angular momentum conservation law it follows, that the impact parameters are 
also equal 

b = b' (38) 

for the collisions; thus, the equality ( 136]) follows directly from that the Jacobian determinant 
of the transformation of not primed velocities into primed velocities equals to unity; the 
equality does not depend on form of the functions 0, /, the only requirement is that integrals 
are determined and converge - cf. with [2J], Chapter 3, 5 5, Section 3. 

ri 

Following Hecke's idea (see Chapter XXII), we multiply both sides of the equation 
( 135]) on In f-, integrate over c?, sum over i and transform integrals, taking ( 121]) . (136]) into 
account; as a result we obtain: 

(AO) AO) \ 
44* (4 o, /i o) -r/r)%*^^=°- < 39 > 
h ■> j ) 

Integrands in (135]) can not be (strictly) less than zero, therefore the sum in (135]) can be equal 
to zero under the condition only, that all integrands for all values of integration variables 
are vanished (all considered functions are assumed to be continuous in each point of their 
definition domain), i.e. 

f{0)i j{0)l ^ y(Q) j(0) ^ 
1 3 % J 



or 



In /f' + In if' - In if - In if = 0. (41) 
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Hence, In should be expressed linearly in terms of summational invariants of collision 

4 1] = m h (42a) 
V? } = rrna, (42b) 
4 3) = \m lC l, (42c) 

for a collision of the i-molecule with the j-molecule conservation of the collision invariant 
ip*P is expressed by equality: 

jf,®' + _ jf,® _ = (/ = 1, 2, 3). (43) 

Linear combination of the summational invariants ib{ is an additive invariant. There 
are no other summational invariants, that depend on velocities of molecules and are linear 
independent of if>i . six scalar unknown quantities (six components of velocities of the 
molecules after the collision ci, ci) are completely determined by six known components of 
velocities of the molecules before the collision c-, c* by two free geometric parameters, that 
define a collision, for example, such as impact parameter b and azimuth angle e (see above), 
that gives us two constraint equations, and from four scalar equations (|4"3"]1 . corresponding to 
the conservation of energy and three components of momentum; Existence of summational 
invariant of collision, depending on velocities of molecules and linearly independent of 
would give similar to ( 143|) redundant constraint equation for velocities of molecules after 
collision and velocities of molecules before collision. 

Therefore, 

ln/ (o) = a (i,o) + a (2 ) o) . m{C{+a <M± m{< * {i (44) 

where a^' ' and a- 3,0 '' are some, independent of c-, scalar functions of spatial coordinates 
(that are defined by the radius vector r) and time t, and a- 2,0 '' is a vector function of r and 
t [as it follows from the equations fj4*3|) for 1 — 2,3, the functions a?'°' ) and a- 3,0 ^ are the 
same for all inner components of mixture]. Or 

In fP = In 4°' 0) + a ( - 3 ' 0) - m- fa + /a (3 ' 0) ) ' , (45) 

where <yf'°^ is a new scalar function of r and t. I.e. the general solution of the equations 
system fl3"5]) can be written as a set of the Maxwell functions: 

3/2 m,(c-0( 2 -°)) 2 

e(n, .fj.n, I mi \ ~ 1 2 k p^ 



f = ri7il - V I'^l, (46) 
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where k is the Boltzmann constant, 

3/2 



(3,0) 



/3 (2 ' 0) = ~F07, (48) 



(2,0) 



cr- 



(3,0) ' 



Particle number density n, of the i-component, mean mass velocity u and temperature 
T~ of inner components of mixture are introduced by definitions: 

rii= / fidci, (50) 



fcT *I>=£ jlrn i (c i -u-) 2 f i dc h (52) 



in [52] A: is the Boltzmann constant. From (150~]) -( I5"2"]) the equality is obtained: 
3 
2 

ieJV ieAf ieTV 

that is convenient to use below instead of definition (|52p . 

According to definitions fl50|) . ([5T]) . fl53|) . in addition to the asymptotic expansion (J2j 

it is necessary to determine asymptotic expansions for particle number density rij of the 

i-component 

^ = nf ] + 0nW + #\ (2) + ■ ■ • , (54) 

mean mass velocity u 

u~=u {0) + 6u { } ] +6 2 u {2) + ■■■ (55) 

and temperature T~ of inner components of mixture 

T- = T- (0) + OT^ + 6 2 T {2) + ■■■ . (56) 

Substituting (1271) and (154|) - (!56|) in (|50|) . (I5T1) . (|53|) and equating terms of the same infinites- 
imal order we obtain Card ( N ) +4 scalar relations, that connect asymptotic expansions (J2j 



and (El- (ESI): 



= <6JV, (57) 
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/r r 
m i c iff )dc i = Y m i ( n i w ") W = E m iJ2 n i~ S)ui ~ S) = E P {r ~ s)ui ~ s \ ( 5 

;_V l£jV i&N s=0 s=0 



E / \ m ^Vt )dc t = \ k Y (^-) (r) + \ E m i( n i u ') ' 

- k Y n (r - s) T^ + -YY P (r ~ s)u[S ~ q) ■ uiq) - (59) 



(r—s) (s—q) (q) 



2 2 

s=0 s=0 q=0 



In (I51?j) the notations are introduced 

p(-) = 5>,nf" s) , (60) 



E4" s) - ( 61 ) 

In particular, for r = from fl57j) - fl59l) we obtain expressions for arbitrary functions 
/3i lfi) (r, t), /3 ( 3 0) (r, t) and /3?' 0) (r, t) in gSJ through the zero order approximations to local 
values of the i-component number density, the mean mass velocity and the temperature of 
inner components of the mixture: 



0(1,0) 

% 


(r,t)=nf (r,f), 


(62a) 


M 2,0) 


(r,t) = u {0) (r,t), 


(62b) 


^(3,0) 


(r,t)=T! 0) (r,t). 


(62c) 



According to fl22l) . fl32l) . zero order integral equations, from which the velocity distribution 
functions fl - 1 of particles of outer components of the mixture are found: 

^(/f } ,/ (0) ) = // (fPf^-f^f^)kidkdc = (ieN), (63) 

- are simpler than equations (|35|) and differ actually from (1351) only by lack of summation 
over components. Therefore, similarly (|46p . the general solution of the equations system 
(|63|) can be written as a set of the Maxwell functions: 

3/2 ^-4 2 '°T 

2ttA;/3^ " 



/f = ^ 0> (r^y) •"^^ (64) 
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where pi 1 ' ^ and /3?' ' ) are some, independent of c-, scalar functions of spatial coordinates, 
defined by the radius vector r, and time t, and /3l 2,0 ' ) is a vector function of r and t. 

Let's add to the definition of the number density of particles of i-component definitions 
of mean velocity u- and temperature T- of outer component of mixture: 

u-n-m^ / rriiC-f-dCi, (65) 



| kTj m = J ^ m- (c - U|) 2 /i rfc- ; (66) 
from (150]). (J55]) . (El]) the equality is obtained: 

~ kT- n- + i m? n • m ■ = ^ i c| /=- dc- , (67) 



that is convenient to use below instead of definition f[66|) . 

Let's enter similar (|55l) -( l56l) asymptotic expansions of outer i-component mean velocity 

ui 

u- = uf + duf + 9 2 uf + ■■■ (68) 
and outer z-component temperature T- 

T- = T, (0) + OTP + 6 2 T? ] + ■■■ . (69) 



Substituting (|27|), (j54jh (|68p, f|69f) in (|50p, (j65]h (j67|) and equating terms of the same in- 
finitesimal order we obtain for each i 5 (scalar) relations, that connect asymptotic expansions 

(EZD, PI, PD, PI): 

'jf<*c< = nf\ (70) 

/r r 
mic-f^dc- = m { (n^) (r) = nl r " s) nl s) = ^ pl r_s) ul s) , (71) 

s=0 s=0 

■ m^ffda, = 3 -k(n { T { )^ + ±fnf(f*i?) (p) 

= 2 A, y nl r - s) T, (s) + -m, v y n { r s) u { r q) . 

s=0 " s=0 q=0 

= § * e » w 4 e e *r"«r* ■ 4". (72) 

s=0 s=0 <?=0 
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cf. with (E7D-fl59D. In (E]), (H2D the notation is used 

p\ = rri{n\ . {73) 

For r = from (170 1) - (1721) we obtain expressions for arbitrary functions fifi (r, t), f3f' ^ (r, t) 
and /3r 3 '°' ) (r, t) in flMl) through the zero order approximations to local values of the number 
density, the mean velocity and the temperature of outer i-component of the mixture: 



(1,0) 
i 


(r } t)=nf< 


[r,t), 


(74a) 


(2,0) 

i 


(r,t) = uf 


(r,t), 


(74b) 


(3,0) 
i 




(r,t). 


(74c) 



pi 

For r > 1 the velocity distribution functions of inner components of gas mixture are 
found from the integral equations system ( 13TT) . which, taking ( |23i) and ( l4"0]) into account, 



can be rewritten in the form 

r— 1 r— 1 



jgTV 5=1 JeTV «=o 

=-EM'W/f)-EM>r.w) 



jeN jeN 



Efff^fPi^ +Xf ) -Xf )/ -Xf)%^dc, (ietf), (75) 



in (175|) functions /r are written as /?"' ) = fl X* > where are new unknown functions. 

The left-hand sides of equations (175|) involves functions, that are known from the previous 
step of the successive approximations method. Unknown functions X? enter, linearly, only 
into the right-hand sides of equations (|75p . Therefore the general solution of the system 
of equations (13T1) is a family of functions of a form {f~ = 5? + £^}i e /V) where {si 7 ^ = 
^1 }j- e jv> {^'"' ) — JieiV; a family of functions {^r }j f 6 # is a particular solution 

of the system of inhomogeneous equations ([75]) and a family of functions {</>^}i e ft is the 
general solution of the system of homogeneous equations 

= E // fffT Of + 4 ] - 4 V - % rfk rfc ; *) • ( 76 ) 

Multiplying equations (j76p by </>/", integrating over all values of c-, summing over i and 
transforming integrals, as it has been made in deriving (1391) . we obtain 

jE jjlflf > {*f + *f-*T'-*f)\^ l lc t ic 1 =a. (77) 
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From (I77|) we conclude, cf. with ( 139]) and ( H4|) . that <fr? are linear combinations of the 
summational invariants of the collision ipf* (I = 1, 2, 3): 

j» (1,0 i (2,r) . (3,r) 1 2 />7 \ 

where a I 1 and a- 3 ' are some, independent of c 4 -, scalar functions of spatial coordinates, 
defined by the radius vector r, and time £, and a- 2,r ^ is a vector function of r and t (as well 
as above, arbitrary functions ol?'^ and oA 3 '^ are identical for all inner components of the 
mixture), and, hence, 

tf> = /f(a?" + aW. W+ a^^) (i 6 iv). (79) 
To simplify further evaluations according to the expression for f~°\ see (pE5)) and (152)1 . let us 



rewrite ( 179)) as 
# = if 



+ /3 ( 3 r) • m { (c; -vf®)+ P {3 > r)1 -m; (c { - 2 



(ieN), (80) 



where /^ 1,r) , /3 (2 ' r) and f3^ r) are new functions of r and t. Family of functions {x^ /ikft i s 
a solution of the system of inhomogeneous equations 

F f ) = E// 4 0) i/ 0) (4° + Xf - - xf) % ^ d C/ (< g n) , (81) 

where F? denote left-hand sides of the equations (1751) . taken with opposite sign. 

Multiplying equations (J5TJ) by ip? (I = 1, 2, 3), integrating over all values of c- and trans- 
forming integrals as above, we obtain, taking ( H3)) into account, as necessary condition for the 
existence of solutions of the system of integral equations ( 1811) . the necessity of the fulfillment 
of equalities: 

|vfifW=0 (<e#), (82a) 
^/vfifW^O (Z = 2,3). (82b) 

Among (infinitesimal) set of particular solutions of the system of equations ( 18 ip . different 
from each other on some solution of the system of homogeneous equations ( 175)) . unique 
solution {#i j-j-gjv ma y ^ e chosen such that 

^fi 0) ^dc{=0 (ieN), (83a) 
E/4°4 0) ^W=0 (1 = 2,3). (83b) 
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Having substituted expression for ff ] [teN 



J i 1— 'i 

■'i i ■'i 



o(l,r) , a (2,r) / (0)\ . o(3,r) 1 / (0) x 2 



34) 



in (l5Tj) - (159"]) . taking (SBJ, (IS0|) - (152|1 and ([S3]) into account, we obtain a system of Card ( N) +4 
algebraic equations [constraint equations for asymptotic expansions ff27|) and (jS"4"|) -( l5"6"|) ]: 



i e N 



(85) 



u (0) £ m-nf^ r) + p^kT^(3^ r) + - p^kT^u^^ r) = £ (86) 



n i0) 3kTi 0) + mAu {0) 
i 1 v 

i'eJV 
3 



4 



(0) 



ft 



(3,r) 



s=0 



(87) 



s=0 g=0 

from which we find expressions for functions /?i ( r ,t) } /3? ,T ^ and /3~ 3,r ^ (r*,t) through 
(variable) coefficients of asymptotic expansions of the particle number density of i-component, 
of the mean mass velocity and of the temperature of inner components of the mixture 



(i,0 



n\ 3 l 



n (o) 2 n( o) T (o) 
1 1 



^ ^(r-.) r (-)^ _ ft (r) T (0) 



2 n(°)£;T~ (0) 



s=0 

r s 



'it- ■ u"' — p y 'hi 



s=0 g=0 



1 



n(°)£;T~ (0) 



it 



(o) 



,s=0 



{r-s)„As)\ _ -(r),.(0) 



(88) 



/3 



(2,0 



p(°)fcT. 



(0) 



s=0 



(89) 
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(3,r) 



k 



s=0 



1 1 

+ 3 — 



(n^T-W) - ftW# 



n(°) [ fcT~ (0) 



p v "Ui • U- ' — p y Mil 



_s=0 q=0 



3 n(°) r fcT. (0)N " 



it 



(o) 



s=0 



p v "IT- I — p y U 



(90) 



Then the fulfillment of equalities (182]) can be considered as the differential equations, the 
r-order equations of gas dynamics, for finding n{ 1 \ u{ , T^ r ^ (r = 1, 2 . . .). 
For r = 1 from (l88D-(l90l) we have 







3 T. 



(i) 







(2,1) 



nf 2 T (o) 



fcT. 



(o) 



(3,1) 



1 T. 



(i) 



(91) 
(92) 
(93) 



The partial solution of the system of inhomogeneous equations ( 1HT1) {$* satisfying 
may be constructed, for example, using expansion of $t (c?) in series in terms of 
Sonine polynomials with expansion coefficients, depending on r and t (see (2] or 10]); such 
construction proves existence of solutions of the s yste m of integral equations (I75|) (one could 
simply use here and above Fredholm's theorems (l^ . [l^|). 

(r) 

For r > 1 the velocity distribution functions of outer components of gas mixture p are 
found from the integral equations system (1331) . which, taking account of (1221) and analogous 
PIP equality 

y(0)/ y(0)/ = j(0) y(0) 



(94) 



can be rewritten in the form 

r-1 



" r 1 + E * (ft*. /"») + {ft 1 -", it 

s=l j-Li s=0 

= -J i {f^\f (0) )-J i {ft\f m ^ 

/fV (0) + X {r) - - kfdkdc (iEN), 



(95) 
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in f )95|) functions fr = p X» , where are new unknown functions. 

Equations (1951) differ from equations (!75|) only in the left-hand sides, that are known from 
the previous step of the successive approximations method, and in the absence of summation 
over components in the right-hand sides of equations (1951) . Therefore similarly to, how it 
has been done above for inner components, we obtain expression for the general solution of 
the system of homogeneous integral equations, corresponding to (l95j) , 



dr) (1,0 , (2,r) . (3,r) 1 2 



Hence, 



s-(r) p(0) I (1,0 i (.-./■) i.-). /•)■■- ■> 



.,(2.0 



(3,r) 



1 



(i G N) 



(96) 



(97) 



or 



(IeN). (98) 



Necessary condition for the existence of solutions of the system of integral equations (f9"5]) 
can be written as 



J 1>f>FPdcr=0 (ieN, 1 = 1,2,3), 



(99) 



where Fr denote left-hand sides of the equations (1931) . taken with opposite sign. 

Among (infinitesimal) set of particular solutions of the system of equations ( |95|) by the 
condition 

/ ipfff ] ^pdc; =0 (i G N, I = 1, 2, 3) (100) 
unique solution may be chosen. Having substituted expression for (% G N) 



f (0 = s (r) , t(r) 

= f (0) «pW , f (0) 



$ 1,r) + M 2 ' r) ■ ( Cf - «r ) + /T"^ Met- ^ 



(0)\ . fl (3,r) 



.(0) 



101) 



in (1701) - (1721 . taking ( 173"j) - (171|) and (11001) into account, we obtain for each index i a 

system of 5 algebraic equations [constraint equations for asymptotic expansions (127]) and 
(El, (El- (Ell: 



(102) 
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« 1 l 1 l 'l 1 ' 2 2 4 1 ^ 1 

s=0 



(103) 



-n 



(o) 



3*T, W +W« (0) 







(l,r) 



i i % i 



+ -kT^ 

4 4 



5n (0) fcT, (0) +p (0) f M l 0) 

2 2 ' J V I 



/9 



(3,r) 



(r— s) (s— g) (g) 

Oi Ui ' ■ Ur' 



(104) 



s=0 



s=0 g=0 

from which we find expressions for functions /3i ( r ,t) } ftf'^ (f,t) and f3i 3 ' r ^ (r,t) through 
(variable) coefficients of asymptotic expansions of the particle number density of i-component, 
of the mean velocity and of the temperature of i-component 







(r) 

n) ' 3 1 



n (0) 2 n (0) T (°) 



s=0 



Ht(»)1 



2 « 



2 n (0) fcT, (0) 



/ -r / -r ' J 2 2 2 \ % 



s=0 q=0 
r 



+ 



,(0)^(0) u i 

i i 



(0) 



Or Ui 



s=0 



(105) 



(V) 



p (0) fcTr (0) 



e (pf- ) -- w 



M (o) 

Ur | — Pr Ui 



s=0 



(106) 



(3,r) 



+ 



3 ni 0) (kT^ 2 



s=0 



/ ■* / -i i i i r i V i 



s=0 g=0 



3 n (o) A T (o)-' 



U 



(0) 



£(p^uf)-pfuf 



(107) 



The fulfillment of equalities fl99l) can be considered as the differential equations, the r-order 
equations of gas dynamics, for finding ni r , itl r-1 \ Tr T ~^ (r = 1, 2 . . .). 
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For r = 1 from fl 1 5 [) - fll 71) we have, in particular: 



1 i 



Pi 
1 i 



01 
1 % 



(3,1) 



(1,1) 



(2,1) 




(109) 



(110) 



(108) 



kT (o) T (o) ' 



i i 



cf. with (EI])- 093]). 

The partial solution {#r °f the system of inhomogeneous integral equations (1931) . 

satisfying (llOOp . may be constructed similarly to, how satisfying ( 1831) partial solution 
lieN °f ^ e s yst em °f inhomogeneous integral equations ( IHTl) is constructed, see above. 

IV. THE SYSTEM OF FIRST ORDER EQUATIONS OF MULTICOMPONENT 
NON-EQUILIBRIUM GAS DYNAMICS 

Let us consider in more detail the system of infinitesimal first order equations ( l82|) . ( 1931) 
(r = 1), derived above as the necessary (and sufficient) condition of the solution existence 
of the first order integral equations system ( 1751) . ( 1931) (r = 1). 

To simplify transformations, according to the expressions for velocity distribution func- 
tions of particles of infinitesimal zero order (H3j) . ( )63|) . functions , may be used in 
( )82|) and ( )99|) (r = 1) rather than functions ipi , tpr, respectively: 



for inner components C- = c- - u. , for outer components C- = c- — ui . 

At transformation of differential parts of the equations (I8"2l and (1931) we use equalities: 




(111b) 



(111a) 



(111c) 




9* 



(112) 
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yWc. . ^lH dc = — ■ I ¥ l) c- f (0) dc -la- ^ f (0) dc- 

-(0) 



9 — m - (°> d^ (l) 

-w-qT > ( 113 ) 



(l) Xi ,djf d = _ f (_d_ tfi)Xi\ AO), 
1 d Cl dCl J [dc, mJ U dCl 



~a ^ (0) 

= -n~-*f ) — ■ (H4) 

OCi ITLi 

In (jll2l) - (IH4p the bar above symbol with index denotes the average of the value: 

F (o) = l [ Vf (o) dCi . (115) 



rii „ 

r and Cj are considered as independent variables; At averaging in (I114p it is assumed, that 
external force X^ acting on the particle of species i, is independent of the particle velocity, 
it is assumed also, that integrals, depending on external forces Xi, are convergent, and 
product typXifj® tends to zero, when Cj tends to infinity. 

After simple transformations from f )82|) and ( 199]) (r = 1) we obtain following system of 
infinitesimal first order equations of multicomponent non- equilibrium gas dynamics: 

^ • n<°>tt<°> (iEN), (116) 



dt dr 

Q ,„,, ^(n\ f \ /■«■) (9 



p (0) ^- + f • P (0) + V j ( 1 = Vn (0) x,-p(°W 0) .^ w (0) , (ii7) 

dt dr f— ' *>>y f-f * * r 5r 



m + . (0): s^ + E a.^m (118) 



dn {0) Ft 

= • nP'tii * (* r G#), (119) 



¥ + |:-pf + E^ = "f^-4V-f-|:«f (^*), (120) 



dt dr 1 1 dr ^— J <9r 
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In accordance with the general definition of pressure tensor of i-component of gas mixture 

p.S/m.te -«,)(«, -«,)/,&, (122) 
and with the general definition of i-component heat flux vector 

Qi= j \ m i( C i ~ U if ( C * ~~ U i)fi dC i ( 123 ) 

(cf. with [2], Chapter 2, §§ 3, 4) in (fTTR|> -tfI2Tl) 



—/ ^rs \-(°) 

p<°> = nfm t (c- - u (0) ) (c* - u (0) ) = # kT^U = p^XJ (124) 

is inner components pressure tensor of zero order, p^ is inner components hydrostatic 
pressure of zero order, U is the unit tensor, double product of two second rank tensors w and 
w' ([2], Chapter 1, § 3) is the scalar w : w' = ^ Q w a pw'p a = w' : w, 

(0) 



is inner components heat flux vector of zero order, 

1 ~, 7^ (0) q 

£(o) = 1 J- n (o) m| _ u (o) j = 3 ^(0)^(0) (126) 

is zero order internal energy of particles of inner components per unit volume, which is equal, 
in this case, to energy of their translational chaotic motion, however, the energy transfer 
equations, written in form (IllSp and (112 ip can be used in more general well (cf. 

with |l0|, Chapter 7, § 6), in (I124l) - (ll26p averaging (11151) is performed with Maxwell function 
ff from ([46]); 

(o) 



p^nfWk-^ K-ul 0) = nl°te (0) U = M 0) U (127) 

is i-component pressure tensor of zero order, pf* is z-component hydrostatic pressure of zero 
order, 

"(0) 

nf^mi ( cr — ul ^ I I cr — I 



q f = l^ mi (c i -^) 2 (c i -^) =0 (128) 



is i-component heat flux vector of zero order, 



r(0) 



= In^mAci-u^) 2 = |n (0) ^ {0) (129) 
% 2* V 1 / 2 4 * 
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is zero order internal energy of particles of i-component per unit volume, in (I127p -( TT2"9]) 
averaging f l 1 X 5 j) is performed with Maxwell function fr from 

Analytic expressions for integrals Jp°jj, J em from ( I117p . f lll8[) and (I120p . (112 ip are given 
below - see ([£24]) . (1A"28]l and (|B2]), f lB3|) . 



In particular, if mean velocities and temperatures coincide for all component of mixture 
(the set of outer components is empty), then from (I116p - (lll8p we obtain the system of 
infinitesimal first order gas-dynamics equations of the Enskog- Chapman theory [it may be 
noted in passing, that the sums over all indexes i of collision integrals in ( 19 9 p are equal 
to zero in all orders of the method of successive approximations, it is easy to check, using 
equality ( 136]) . from the physical point of view this statement is reduced to, that total impulse 
and (kinetic) energy of particles of gas do not change in their collisions between themselves] : 

<J ' r niV 0) Uefi), (130) 



dt dr 

*° "ir^r- = E «' 0,X -' - ^ ■ Tr "™ (131 > 

+ (132 ) 

at or or Or 

It follows from (T130D . (ri32|) . taking (Tl2l)- (ri26|) into account, that the flow of gas, de- 
scribed by the system of first order gas-dynamics equations of the Enskog- Chapman theory, 
is adiabatic: 



D 
Dt 



aw Mo)V 3/2 



0, (133) 



in the adiabatic equation (1 1 3 3 [) 



V. VELOCITY DISTRIBUTION FUNCTIONS OF INFINITESIMAL FIRST OR- 
DER 

Using the system of infinitesimal first order gas-dynamics equations it is possible to specify 
the form of velocity distribution functions of the first order. In particular, we shall see, 
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that the velocity distribution functions, obtained by the proposed method within Enskog's 
approach, i.e. when outer components are absent, and by the Enskog's method are equivalent 
up to the infinitesimal first order terms of the asymptotic expansion (inclusive). 

As 



In/ 



(0) 



const. + Inn! "* lnT- ^ — 



m- \ c--u 



(0) 



(135) 



see (H6T) and ( |62i) . the system of integral equations (J75i) (r = 1), from which {x?}ieN are 
found, excluding time derivatives by equations f l 1 1 6 [) - nTT8|) . taking f ll24p - fll26p into account, 
can be written in the form: 



jeN 



Of 



dr m„- dc; 



ft 
fP 



dt 



dr 

«91nT W 
dr 



+ 



mt dct 



)+£%(/f./f) 



+ 2(C i C t ):-uV + - M C r d i +s i 



E//4 0> /f ( 



(0) f (0) ( (1) (1) _ (1)/ _ (1)/ 
Kg * Aj X? Xj 



X { P')k^dkdc^ (ieN). (136) 



In (11361) . as well as in the previous section, C- = c- — u 



(0) 



1/2 



m; 



2kT 



(0) 



(137) 



C- is the modulus of vector C «• ; for arbitrary second rank tensor w 

w = w — — U (U : w) 



is tensor with zero trace; 

(0)\ / (0) (0) \ „, m / (0) 



1 dr / \ n(°) 



P 



(o) 



<9r 



p(0)p(0) 



p(o) 



(138) 



(139) 



"p(0)p(0) 



c .- E/S-^(^!-i) E/. 



(o) 
£,0 



Lvj..ff(») f ( 

(o) «i \ J i ,J j 



leNjeN 

(0) ,(0) 



leNJeN 



jeN 



(140) 
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The general solution of the system of homogeneous integral equations, corresponding to 
the system of equations (I136p . was constructed above in the section [TTT1 A partial solution 
lieN °f ^ ne system of linear integral equations (I136p can be sought in the form: 

*?' = ~ A f ■ - b, : |y o) - ««»» £ « • *, - % d«) 

jeN 

Substituting ( 1141ft in (1136j) . we find, that family of vector functions (cf)}* 6 # must satisfy 
the system of integral equations 

/r(<?-§)c,= E // /f'/f (A'+^-AJ.-A^gdkdc, (iEJv), (142) 
family of tensor functions {B-(c-)}- £] ^ must satisfy the system of integral equations 

2/f (C,C ; ) = E // 4 0) // 0) • B j B ;) % ' /k A>) • (143) 

Double family of vector functions {Di (c*)}* ? gJ y must satisfy the system of integral equations 

The- ? * «^ «^ 

2 J,h£N 

It follows from the definition d^, see ( 11391) . that 

E<%=0, (145) 

iev 

therefore for every i one of coefficients ZH may be considered as equal to zero; for symmetry 
of expressions it is usually assumed, that D\ = 0. Using (I145p . the system of equations 
(11441) may be rewritten in the form: 
1 



/fcV(%-%) 



(o) 
n) 



= Y f {0) f { ? ] ( + - Dh - D' - D'f - D'f + D'f + D'f) k;;dkdc; 
J J 13 \ 1 3 1 3 1 3 1 3 J 3 3 

U,h,keN\ (146) 

n 

cf. with [10|, Chapter 7, § 3, (3.32). In ( I146P 5^ is Kronecker's delta symbol: 

,0 (i^k), 

S lk = { K ' (147) 
1 (i = k). 
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Family of scalar functions {S- (cA}- eI ^ must satisfy the system of integral equations 

fj\ = E // ff^T {St+Sj-Si- Sj) % r/k da, (/-:.V). (148) 



Solving of the system of integral equations (11481) (this question will be considered, possibly, 
in our next article) becomes a little simpler, if functions S- are written in the form: 

s i = s r E J l% + S h ( 149 ) 

Conditions of solubility of the system of equations ([132]) . ([Tig]) . (fT44D . ([Tig]) , following 
from (|82[) [in ([82]) must be replaced, respectively, by left-hand sides of equations (11421) . 
(TT43D . ([HID and (TUB"]) ], are fulfilled, because conditions (J82D in the form ([TTB]) -(nT8l were 
used in the derivation of the integral equations system (I136p . 

Because r and t do not enter explicitly in equations (11421) . (I143|) and (11441) . enters 
only in combination (c? — it-°^ = C- (after replacement of integration over c- by integration 
over C ■) and the left-hand sides of equations (11431) are symmetric tensors with zero trace, 
functions {A-ic-)}.^, {B f (c^)} fe # and {Di (cA} iJe $ can be sought in the form: 

A^A^nfA^C?, (150) 

B { =B { (nf,C t ^){C{C { ), (151) 

/>; nU^.c.r^c, (152) 



where A- (nf\ Q, T-°^J , B-yn^\ T^J and Dj (nf\ j are scalar functions of 

n (0) ,QandT~ (0) . 

Imposing on the functions { ^ ne condition (185]) . from (184]) and (]9~T]) -( l9l?]) we obtain 

the following expression for the infinitesimal first order velocity distribution functions of 
inner components particles of gas mixture (i E Nj : 

if - if ( ' ■ " B, = - A« 

+ /f (^)"2^ + ^- + ^) ^) 2 • (153) 

Similarly, the system of integral equations ( ]95]) (r = 1), from which {x^ }iejv are found, 
excluding time derivatives by equations (I119I) - (I121I) taking (I127p - (I129I) into account, can be 
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written in the form: 



dt 



a AO) o f (0) „ o f (0) 

°h , _ i A i gZi i \^ 7 f f (o) f (o) N 

(0) , dln/f dln/f X { din if 



f _(o) 



9i 



<9r 



+ 



mr dcz 



„ 5\ r> In / n o d 

/fV (0) (xf ) + x w - xf } ' - Aydkdc (i € iv) : 



(154) 



where C- = c- — u 



(o) 



1/2 



2kTi 



(o) 



(155) 



£- is the modulus of vector C? 



L r ..v / (o) - — (-c 2 - \\ v / {0) 

(0)°» Z^ J P ,ii (0) I 3 S L JZ^ J E,i 

Pi jjki Pi j^i 



(156) 



A partial solution {#» x }{ e ^ of the system of equations (11541) can be sought in the form: 



(i) 



dhiT$ 0) d (o) 



(157) 



Substituting (1 157ft in the system of linear integral equations (11541) . we find, that family of 
vector functions {A{(c{)}- £ ^- must satisfy the system of integral equations 

/<°> (<?-0C,= JJfPf^(A i+ A-A i -A)k i( n,dc (ieti), (158) 

family of tensor functions {B- (c^)}-^ must satisfy the system of integral equations 

2ff{C l C l ) = If /fV (0) (B/ + B-B;,-B')^dkdc (ieiV); (159) 



family of scalar functions {^(^lieJV mus t satisfy the system of integral equations 



fjPs { = II fPf^^+S-Si-S^kdkdc (ieN). 



c(0) n (0) 



(160) 
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solving of the system of integral equations (I160P becomes a little simpler, if functions S{ are 
written in the form: 

• s '<- s " ■ Y, J Tu • s h ( 161 ) 

Conditions of solubility of the systems of equations (I158p . ( I159p . (I160p . following from 
(199|) [in (J99J) must be replaced, by left-hand sides of equations fTT58|) . (11591) and (TT60]) ]. 
are fulfilled, because conditions ( 1991) in the form ( 1 1 1 9f) - ( TT2TT) were used in the derivation of 
the integral equations system (11541) . 

As r and t do not enter explicitly in equations (I158p and ( I159p . uf^ enters only in 
combination ycj— u^j = C- { (after replacement of integration over c by integration over 
C) and the left-hand sides of equations (I159P are symmetric tensors with zero trace, functions 
{A - (c*)}- e fi and {B-(ci)}- e] ^ can be sought in the form: 

Ai = Ai(n?>,C{,Tf>)c { , (162) 



Br^Br{vf,C h Tf){C^Ci) (163) 



where A, ( n^, Q, ) and B- \ are scalar functions of nf \ Q and T? ^ 



Imposing on the functions {$1 the condition ( llOOp . from (llOip and (ll08p -f fTTUj) we 

obtain the expression for the infinitesimal first order velocity distribution functions of outer 
components particles of gas mixture /? (i G N) : 

\ dr dr 

+/f ' b?-2^ + ^' m -' C{ + 2 mfC? ' <164) 

\ i i i i i / 

Let us consider a case when outer components in a mixture are absent. Up to the 
infinitesimal first order terms (inclusive - see {(]], Chapter V, § 2, Section 2, definition 2) the 
expression, see (I153p . 

ff> + /f (^J " ^ + ^ ■ m { C f + - mf C|j (165) 

coincides with asymptotic expansion of the solution in the Enskog- Chapman theory 

3/2 m-fcr-u'. 1 ^ 



^-""'i^' C ^ ' <166) 
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where nl = n? +n? , = and T-! 1 ^ = T- +1^, cf. with the Taylor expansion 

of at the point vfP\ T^j . This statement can be written in the form: 



tr ~ sr +fn±-lw + *®>- miCi + w*w 2 m -' c| 1 • (167) 



Equations (I142p . ( I143j) and ( 1146 j) differ from analogous equations [l0(, Chapter 7, (3.34), 
(3.33) and (3.32) in the Enskog-Chapman theory essentially only in use and T- ^ 

instead of n,, u and T (i.e. n^, lA 1 ^ and T^). Because functions {A^C')}^, {B^c^)}^^ 
and ( c i)}iJeJV are °f the infinitesimal first order, in the first, second and third terms 
in the right-hand side of (1141 jl and also in ( I150p -( ll~52l) functions n|-°', u~ and T- ^ can be, 
respectively, replaced by functions nl 1 ^, and T~ . Therefore up to first infinitesimal order 
terms the expression, cf. with f |153j) . 



coincides with the expression for the function f- 1 ^ in the Enskog-Chapman theory 



Consequently, up to first infinitesimal order terms the functions /I = /? + /? coincide 



with the functions = fi°^ + f? , obtained in the Enskog-Chapman theory: 



/f if- (170) 

As a result, with the same precision (cf. with Chapter 8, § 4 and flo| . Chapter 7, § 4) the 
expressions for the diffusion velocities of components of gas mixture, for the i-component 
heat flux vector, . . . coincide, in particular, the expressions for i-component pressure tensor 
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coincide: 



J 1 ] 



[i] 



1>; ' / '», (c ; - u - ) (c f -ul 1] ) ffdc. 
1. ( n (?hT® + n^kT^ + nfkT^) U 

(0)\ / (0)\ /.(O) D 

* 1 ' c/ - u\>) f 'B. 



m- c- — u 



9 (o) 



^nM 1] U-2 



1 mf 



15 2A;Tl 



4 <9r 1 



[i] 



d_ 
Or 



u 



[i] 



(171) 



cf. with incorrect opposite assertions (about Hilbert's method for the kinetic Boltzmann 
equation asymptotic solution, which, in principle, should lead to the same results, as the 

3, B, 3; in (QHD the 



proposed in this article method, see section |IT]), for example, in 
notations are used: fi- is viscosity coefficient for i-component, for arbitrary second rank 
tensor w 



1 



>a(3 2 

is symmetric tensor, corresponding to it. 



(172) 



VI. TURBULENCE AS MULTICOMPONENT GAS DYNAMICS 

As is well known, laminar flow becomes turbulent flow, when some parameter character- 
izing the flow, namely, Reynolds number 

Re = EH^L > l. (173) 

In (I173p . p is the density of gas, u and L are some characteristic macroscopic velocity and 
linear size of the flow, p is the coefficient of viscosity. Having rewritten (I173P as 

Re= P -^, (174) 

cf. with the expression for viscosity tensor in (11711) . the Reynolds number can be treated as 
the ratio of the macroscopic momentum flux, proportional to f^°\ to the viscosity- induced 
microscopic momentum flux, proportional to (from the point of view of the kinetic 
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theory of gases instead of Reynolds number it would be more natural to consider the ratio 
Re = f(°> / fW). Roughly speaking, viscosity, "aligning" the gas molecules according to 
a Maxwellian distribution with the same mean velocities and temperatures for different 
components of mixture, can "process" only the microscopic momentum flux. But if the 
macroscopic flux exceeds the microscopic, the gas flow, necessarily, has begun to stratify on 
components [and other mechanism, leading to an equilibrium state of gas, "turns on" - see 
integral terms in (I116l) - fll2ip . however, from the point of view of the kinetic theory of gases 
it is the same mechanism, collisions of gas particles with each other]. The flow stratification 
on components can be also caused by external factors. 

The gas-dynamic equations system of the second approximation order (i.e. with diffusion, 
viscosity and heat conduction) of the Enskog- Chapman theory, in principle, can not describe 
gas flow with developed turbulence, as 

f(0) 

Re ~ i > oo, (175) 

/(i) ' v 1 

when /W — » 0; with increasing Reynolds number, according to this system of equations, 
the entropy production decreases until the complete cessation of growth of the entropy 
in the gas, see the adiabatic equation (I133p (terms of the gas-dynamic equations system 
of the second approximation order of the Enskog- Chapman theory, additional to the gas- 
dynamic equations system of the first approximation order of the Enskog-Chapman theory, 
corresponding to the transition of gas to an equilibrium state and, hence, to the increasing 
of its entropy, are proportional to f^ 1 '), whereas the entropy of gas grows with growth 
of Reynolds number in experiments, i.e. theoretical (Enskog-Chapman) and experimental 
dependencies of entropy of gas on a Reynolds number are different. 

Mixing is absent in the gas-dynamic equations system of the first approximation order of 
the Enskog-Chapman theory f ll30[) - fll32j) . because all components of gas mixture, according 
to the mass transfer equations f ll30p . move with the same mean mass velocity of mixture. 
The substance from some physical region (even with moving boundary) cannot cross region 
boundary in any way. Mixing in the Enskog-Chapman system of gas-dynamic equations, i.e. 
a possibility to gas particles to cross boundary of physical region, as well as other mechanisms 
(viscosity and heat conduction), leading gas to an equilibrium, and, hence, increasing entropy 
of gas, appears only in the following infinitesimal order of asymptotic expansion of velocity 
distribution functions of gas particles and is related to diffusion velocities (~ f^) in mass 
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transfer equations. The ratio of diffusion velocity of some component of gas mixture to mean 
velocity of this component ~ 1 /Re ~ fj^ j fj® tends to zero, when — >■ 0. Therefore it 
is also not possible within Enskog-Chapman theory to describe observed intensive turbulent 
mixing in gas-dynamic flows with great Reynolds numbers. 

If gas-dynamic equations do not describe turbulent gas flows, then either something has 
been missed during the transition from the exact solution of the system of kinetic Boltzmann 
equations to its approximate solution (by the Enskog method) and next to the gas-dynamic 
equations, or the system of kinetic Boltzmann equations does not describe turbulent gas flows 
and requires replacement. However the last, i.e. necessity of replacement of the system of 
kinetic Boltzmann equations on another system of kinetic equations in transition from gas 
laminar flow to turbulent gas flow, seems ill-founded. 

The gas dynamics of the components with the velocity distribution functions, close to 
the Maxwell functions, with different mean velocities and temperatures, should be described 
by equations (j!16p - (jl2ip . From this point of view, the observed chaotic character of the 
turbulent flow is similar to the chaotic character of the Brownian motion. They differ in scale: 
in the Brownian motion that particle moves stochastically, whose mass is comparable to the 
mass of other gas molecules, whereas in the turbulent flow that body moves stochastically, 
whose mass is comparable to the mass of separate gas components. In flll7l) - flll8l) . (I120p - 
(11211) the integral terms (proportional to n i: rij) can be huge, it explains unexpected (for 
those, who tries to describe turbulent flow by the gas-dynamic equations system of the 
second approximation order of the Enskog-Chapman theory) power of turbulent effects. 



Appendix A: Calculation of collision integrals 

In this section we are dealing with calculation of definite multidimensional integrals 

(if'if' - if Vf) 9ij b db de d Ci dc r (Al) 

,(1) _ ™ ,T,( 2 ) _ ™ r> ,Tr(3) _ 1 ™ n2 




In dAU) = rm, *r = mCi, = 3 miCf, C t = c t - u t ; 



is the Maxwell velocity distribution function of i-component particles, the prime in the 
distribution function implies, that the distribution of particles velocities c' after the collision 



rrii 



3 /2 mijci-uj) 2 
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is considered; to diminish a little inconvenience of notation, the upper index "(°)" at n.j, Ui, 
Ti is omitted. The other notation is specified above. 

According to (|36|) . the integral (lAlj) can be transformed as follows: 




*i° (fTf^-fPf^tobdbdedvdc, 





(Q/(o) Ao), 




vrfrfrgijbdbdedddcj 

- *f) ./:"' /;;"//,; & rf& rfe dc, dc,, (A3) 

As the particle mass is conserved in the collision, integral ( 1A3j) vanishes for \E'- 1 ^ = m^. 
In two other cases, generally speaking, this is not true. 

Statements of the two following simple propositions are used below several times. The 



proposition [2] is taken from [6J], Chapter II, § 1, Section 5. Regulated functions (see 
Chapter II, § 1, Section 3) in formulations of propositions can be replaced with more known 
continuous functions (any continuous function on R is regulated). As complete normed 
vector spaces over the field of real numbers R further in the article vector spaces R or R 3 
over R are considered, with the usual modulus of a real number or a three-dimensional 
vector as norm. 

Proposition 1. / is assumed to be a regulated function on R with values in R, w e R 3 be 
a fixed nonzero vector, n e R 3 be a unit vector. In this case 

7T 

*2ttu) f 

f(wn)n dtt n = f(w cos (9)) cos (9) sin (9) d9. (A4) 

n n o 

In the left-hand side of the integral is taken over all directions of vector n, w ■ n is the 

scalar product of vectors w and n. 

Remark. If w is the zero vector, then the right-hand side of ( L4^| ) may be set equal to 0. 

Proof. Select the system of spherical coordinates, such that the polar axis direction be the 
same as the direction of the vector w. Resolve the vector n into two components: parallel 
n\\ and perpendicular n± to the vector w, - 

(w ■ n) w , , „ 

n = ny + n± = f V n±. (A5) 



w 2 
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Having substituted expression (1A5j) for the vector n into the left-hand side of flA4j) and 
integrating over the azimuth angle, we obtain the required equality, as in the integration 
over the azimuth angle the term, containing n±, vanishes. □ 



Proposition 2. E and F is assumed to be two complete normed spaces over field R, u be 
a continuous linear map of E into F. In this case, if f is a regulated function on interval 
I G R with values in E, then uo f is the regulated function on I with values in F and 

J u{f{t))dt = u[j f(t)dt\ . (A6) 

a \a / 

Proof. Equality ( 1A6I) follows immediately from the expression for the derivative of composite 



function u o /; details of the proof can be found in 



6], Chapter II, § 1, Section 5. □ 



Major difficulties of evaluation of integral (IA3j) are caused by that parameters of Maxwell 
functions for the i-th and the j'-th components are not equal: 

u, ± Uj , T % ± Tj. (A7) 

As result, it is not easy get rid of scalar products of vectors in the exponent (it is desirable 
to have the expression for the exponent as simple as possible). 

As the scatteringangle depends on the modulus of relative velocity of colliding particles 
{see, for example, |2j, Chapter 3, § 4, Section 2 or [10], Chapter 1, (5.26)}, it is natural 
to proceed in (IA3I) to new variables - center-of-mass velocity Gij and relative velocity of 
colliding particles g^, which are related to particles velocities c« and Cj by 

c i = G ii + — -j—gy, (A8) 

TTl ■ 

Cj = Gij f— o -> ( A9 ) 

777,j + m i 



,2 



- cf. with [2], Chapter 9, § 2. For further simplification of the exponent vector Gy can 
be replaced by vector Gy resulting from G^ in an arbitrary affine transformation, which is 
a composition of shift, homothety (multiplication by a scalar) and rotation. The rotation 
arbitrariness is reduced to the freedom in choosing of direction of the polar axis in the 
transition to the spherical coordinate system. Similarly, the vector can be replaced 
by the vector g^, resulting from g^ in composition of arbitrary homothety and arbitrary 
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rotation. The shift of the origin of the vector o 2? - would lead to parametric dependence of the 
final integral on vectors Ui and Uj (cf. with [16], Chapter 3), that is undesirable, as integral 



3]) is supposed to be reduced to Chapman-Cowling integral Flf-^ [see |2|, Chapter 9, § 3 



i.i 



(3.29) and [lOj, Chapter 7, (4.34)], depending on the modulus of the relative velocity of 
colliding particles only. 

In view of the aforesaid, make the following substitution of variables Gij and gr,-: 

9 lj = zig ij , (A10) 

. — w ■ H - w ■ 

Ga = z 2 G tj + z 3 ~ 9lJ + -^g-^. (All) 

In (lAlOfl - flAllj) the scalar factors Z\, z 2 , and z 3 are selected from the condition that the 

~2 ~ 2 

coefficients of and G^ in the exponent be equal to 1 and the coefficient of the scalar 
product <7jj- • Gtj be equal to (cf. with the method of separation of variables): 



TfliTflj 



2(T t -Tj) r ™<™J (A14) 

mj + mj y 2 (m^Tj + rrijTi) 
Analogous substitutions of variables can be used in more complicated situations, for example, 
discussed in [l6J|, Chapter 3. 

With new variables the exponent can be written in the following form: 



g\j + G\j + a w 2 + aiijij ' w + a 2G\? ■ w 



(A15) 



where 



wj = (A16) 

a = -^ + -^, (A17) 




^1 = -^/ T ~~~777 ) (A18) 



"2- (- — ) J t - (A19) 

It is easy to see, that by the above-specified transformations of variables only, without 
using the shift of the origin of vector g^, it is impossible to get rid of the constant term in 
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exponent (1A15j) and, hence, of the constant exponential factor, which will appear hereafter 



in all expressions containing integrals of form flAlj) . flA3j) . Such factors are missing in 



3j, 



Determine Jacobian of transformation of variables (cj,Cj) — >■ (^g^^Gi^j [see flA10[) - 
(jMU)]: 

d(ci,Cj) d^Cj) d(g ij ,G ij ) _ 3 3 <9(cj,c jy 

Z l Z 2 



_ r 3 r 3 ^ ( C » C j) _ r 3 r 3 ^ ( C » C j) _ r 3 r 3 / A nc\\ 

— Z l Z 2 7T7 r — z l z 2 JT( 7 — z l z 2- {AZU 



Now consider the case, when tyf' = NE^ 2 - 1 = rrii (cj — Ui). In view of ( 1A15I) . ( 1A20I) . ( 1A10I) 
(I A 1 1 [) and the equality, following from the definition of k above, 




rm « - ct) = fa - g v ) = -2 fa ■ k) k (A21) 

mi + rrij J J rrii + rrij J 

integral ( 1A3I) can be rewritten as: 

rrii {c[ - f- 0) fj 0) gij b db de dc { dcj 

x exp (- tgfj + G\ + a w 2 + a x g {j ■ w + a 2 Gij ■ tuj) 

x fa b db de dG^ dg {j . (A22) 

Integrating over e in (1A22f) (with fixed g^ and Gij), we resolve vector k into two components: 
the ones parallel and perpendicular to vector g^ - cf . with the proof of Proposition [1] 

J fa ■ k) k de = 2vr cos 2 ( — y^J g tj = ir (1 - cos x) 9iy (A23) 

When integrating over G^ and directions of vector g^, we use Proposition [TJ As a result 
we obtain 



t(0) 

p. y 





m.C, (if 'if ' - if Vf) ^ 6 d& de d Cl dc 3 
rrii (c! - /f // 0y b db de dc { dc 3 

2 



rriiTj + irtjTi w y/n 



-\Qn i n j — V — e "N r i+ m i T < 

(rrij + rrij) w t, 2 

x e~ 9 ^ [giji cosh (g^) - sinh {%£)) ~g% (1 - cos x) b db dg {i . (A24) 



39 



In (1A24|) 



C, = a\w, 



(A25) 



factor <2i is determined by formula (1A18I) . It is easy to check, that the singularity at £ = 0, 



which is possible when w = 0, is actually absent in the right-hand side. Expression (1A24|) 
differs substantially from Struminskii's expression |3J], (8). 



The case, when \l/ 



(0 



(3) _ 1 



rrii (q — Ui) 2 , differs from the just considered one in the 



factor of the exponent in the right-hand side of (1A22I) . Transform difference ty®' — 
according to ( 1A10I) . (I A lip and [2], Chapter 3, (4.9) and taking account of that only the 
relative particles velocity direction changes during the collision (g^ = g'^): 



(3) = m 

i 2 
rrtirrij 

rrii + rrij 
-2zi 3 



[C- - Ui? - (Cj - Ui) 2 



(g'a - 9ij) ■ (Gij - Ui) 



rn 



1 / „' 



C'i - d) ■ (C- + C; - 2Ui) 



(hi • k ) k • { z * Gij + z 3 g 



Ui — Ui 



(A26) 



With respect to its arguments the scalar product is a bilinear continuous function, there- 



fore Proposition [2] can be applied. After integration over e, similarly to (1A23|) . we arrive 
at: 



-2zi 



miiiij 
rrii + m j 

-2n z\ 



(9ij ■ k) k • <z 2 + z 3 g 



Ui — Uj 



de 



rriimj \ ~ ) 7^ ~ Ui-Uj 



mi + m 



3 



1 - cos x) 9ij ■ S z 2 + z 3 - 



Perform the integration over and directions of vector g^, using proposition [Q 



J 



(o) 

E,ij 




~ m,C 2 (/«> 7f - /fjf) 9ij b db de dc t d Cj 




2 



[c' A — Ui) 2 — (d — Ui) 2 



fi fj 9ij b db de dci dcj 



J (0) - u J [ 

u e,ij "i P,iJ 



(0) 



lQ-riiUj e mi T i +m J Ti 



w 



(A27) 



x / / e "tJ { D hij j [<?i,-£cosh (g^) - sinh (g^)} + 2D 2,ij9 2 j sinh (g^) 
x gf L - ( 1 — cos x) b db dg~ij . 



(A28) 
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In (1A28|) : 



D 2 ,ij 



2 rrijTi 
rrii + rrij 



rriimj (Tj - Tj 



2 (mj + m 




(A29) 
(A30) 



The other notations are the same as in (1A24j) . 

It is interesting to note, that for itj = Uj integral (1A24I) and the first term in flA28j) vanish, 
the second term in (1A28j) is proportional (Tj — Tj), that corresponds to energy transfer from 
the "hot" components to the "cold", see the gas-dynamic equations system (lll6[) -( TT2Tj) . In 
view of the sign of a\, flA18|) and definition of £ (1A25I) the first term leads to temperature 
increase when Uj ^ Uj. 



Appendix B: Expressions of collision integrals for interaction potential of rigid 
spheres 



Collision integrals ( 1A24|) and (1A28j) are complicated functions of mean velocities and 
temperatures of separate components, mainly, because of the complicated dependence of 
deflection angle x on (t ne modulus of) relative velocity of colliding particles - cf. with 
Q, Chapter 1, (5.26): 



X(b,9ij) =n -2b 



dr/r 2 



Mr) 



Ij 2 

72 



(Bl) 



In (IBlj) rriij = rriimj /(rrii + Trij) is reduced mass of colliding particles, p (r) - central inter- 
action potential of particles, depending on distance r between them. 

In the simplest case of particles, interacting as rigid spheres with diameters Cj and aj, 
the following analytical expressions for collision integrals have been derived from (IA24I) and 
fTA28l) : 



r (o). 



-riiUj 



rriiTj + rrijTi w ^Jtt 2 
rrii + m j w 2^ 2 13 



-^ 4 2e(^ + 2) + v ^(^ + 4^-4)erf (0 



(B2) 
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J 



(0). 



E,ij 



n inj |p a 2 e"^ 4 [2D 1>ijW ^ + 2) + 2£> 2)ii £ 2 (£ 2 + 10)] 

7T „ 



X 



+ 4£ 2 - 4) + £> 2jij £ (£ 4 + 12^ 2 + 12)] erf (jj 



(B3) 



In (JB2])-(IB3]) o-y = + and notations from flAT6|) . fIXTSli . (jM5l) . (lA29]) - (lA30l) are 

used. 
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